Nucleation Induced Undulative Instability in Thin Films of nCB Liquid Crystals 
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A surface instability is reported in thin nematic films of 
5CB and 8CB, occurring near the nematic-isotropic phase 
transition. Although this instability leads to patterns reminis- 
cent of spinodal dewetting, we show that it is actually based 
on a nucleation mechanism. Its characteristic wavelength does 
not depend markedly on film thickness, but strongly on the 
heating rate. 
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Following several studies on the spreading behavior 
of liquid crystals (LC) from the nCB homologous series 
(4'-n-alkyl-4-cyanobiphenyl) undulative instabili- 

ties have been observed in thin films of the LC 5AB4 
0, and 5CB || near the nematic-isotropic (N-I) phase 
transition. In the case of 5CB, these results have led to 
some discussion whether a spinodal dewetting mechanism 
driven by van der Waals forces is at work, as proposed by 
Vandenbrouck et al. || , or wether the instability is driven 
by a pseudo-Casimir force based on the director fluctua- 
tions in thin nematic films P-|l3[]. In the present paper, 
we show that neither is true for nCB thin films. Instead, 
the instability is caused by textures in the nematic film 
which largely determine the characteristic wavelength of 
the emerging pattern. 

5CB and 8CB were obtained from Merck KGaA 
(Darmstadt, Germany) and Frinton Laboratories Inc. 
(Vineland, NJ) respectively, and used without further 
purification. Silicon wafers (100-oriented, p-(Boron-) 
doped) with a native oxide layer of 2 nm provided by 
Wacker Chemitronics (Burghausen, Germany) were used 
as solid substrates. The wafers were cut to samples ap- 
proximately 1 cm 2 in size and cleaned with a Snowjet™ 
(Tectra, Frankfurt/M, Germany), a cold CO2 stream ef- 
fectively removing particulate and organic contamination 
Jul , followed by ultrasonication in ethanol, acetone, and 
hexane, subsequently. 

Immediately after this cleaning process, LC films were 
spincast onto the samples from hexane solutions. Varia- 
tion of concentration and spinning rate allows to deposit 
films of variable thickness. The preparation procedure 
was performed in a class 100 clean room environment 
at room temperature. Therefore, the films were initially 
in the nematic (5CB) or smectic A (8CB) state, respec- 
tively. Film thicknesses were recorded with an ellipsome- 
ter (Optrel GbR, Berlin, Germany). The samples were 
placed on a heat stage (Linkam THMSG 600, temper- 
ature control better than 0.1 °C) and observed in situ 
with a Zeiss Axiophot microscope equipped with a dig- 



ital camera. Unless otherwise noted, no polarizers were 
used in the microscope setup. 

Observations at room temperature showed films of 
5CB and 8CB to be stable for hours at thicknesses rang- 
ing from 50 nm to 200 nm. Upon heating, a surface 
undulation with a characteristic wavelength can be ob- 
served in both types of samples (see Figure [|) close to, 
but consistently below the N-I transition temperature 
(T NI ). 
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FIG. 1. Surface instability of an 85.9(3) nm thick film of 
8CB at 39.9(1) °C. The scale bar has a length of 100 /Ltm. The 
inset Fourier transform gives a wavelength of 40(20) /im. 

The observed undulation does not lead to a complete 
dewetting of the LC film but rather disappears above 
Tni, such that the film becomes homogeneous again in 
the isotropic phase. Upon cooling, a similar instability 
occurs. Observation under crossed polarizers revealed 
that the undulation is accompanied by lateral 'demixing' 
of nematic and isotropic regions. 

Complete dewetting of the samples occurs only by 
heating a few degrees above Tni, resulting in an array 
of isotropic droplets. As heating and cooling of the sam- 
ple around Tjvj was repeated, as indicated in Fig. |^, 
nearly identical undulation patterns were obtained. We 
found that the instability developed only when the sam- 
ples were heated or cooled at a certain rate. Keeping 
them at a fixed temperature close to the N-I phase tran- 
sition resulted in the nucleation of only a few holes in 
the film. These observations have been made for 5CB 
and for 8CB, at all film thicknesses investigated. The 
pattern was found to vanish faster for thicker films after 
completion of the phase transition (see below) . 
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FIG. 2. Heating/ cooling sequence showing the transient 
instability in a patch of 8CB (scale bar = 25 /im). 



Careful examination of the images revealed the for- 
mation of nematic domains in 8CB films when heated 
above the smectic A-nematic transition temperature for 
the first time after preparation (cf. Fig. ||). This do- 
main pattern is preserved during subsequent heating and 
cooling of the samples which was limited to a few degrees 
around the N-I transition temperature. On the left hand 
side of the figure, contrast is enhanced in order to clearly 
show the domain boundaries, while on the right hand 
side, the undulative instability occurring near the phase 
transition is superimposed on the domain boundary pat- 
tern obtained before on the same spot. It is clearly seen 
from the overlay that the undulative pattern is strongly 
correlated with the domain boundary pattern. 




FIG. 3. Left: 86 nm thick film of 8CB in the nematic 
state. Borders between nematic domains appear as dark lines. 
Right: Picture of the surface undulation superimposed on the 
network of domain borders, emphasizing the coincidence of 
surface undulations during the instability (bright areas, thick 
film) with domain borders in the nematic phase, (scale bar = 
100 fim) 

Since films of 5CB are nematic at room temperature, 
the domain pattern observed immediately after prepa- 



ration is less pronounced, but nevertheless present. As 
in the case of 8CB films, this pattern is not influenced 
during repeated temperature changes of the sample. 

Since we have found that the undulations appear only 
when the sample temperature is swept through the phase 
transition, it is of interest to investigate the impact of 
the heating rate. We have thus varied the heating rate 
from 0.01 K/min to 10 K/min. As shown in Fig. ^, the 
heating rate is indeed a major defining parameter for the 
wavelength of the undulation: The faster the samples 
are heated into the isotropic phase, the smaller is the 
undulation wavelength. 

To explain the instability as such we propose the fol- 
lowing scenario: After preparation (for 5CB) or after 
heating to the nematic phase for the first time (for 8CB), 
a pattern of nematic domains exists in the films that 
is preserved during the course of the experiment, since 
close to the substrate nematic or even smectic order ex- 
ists even at temperatures substantially above the clearing 
point (see e.g. (l5[ ). The domain boundaries act as nucle- 
ation sites for the isotropic phase upon further heating. 
It should be noted that apparently similar observations 
in 5AB4 jjj do not belong to this class of nucleated undu- 
lation phenomena, since in that study the temperature 
was kept constant, and the undulation was not transient, 
but remained once it was formed. 

It is interesting to note that there is no systematic de- 
pendence of the patterns observed with nCB on the film 
thickness. Samples of different thickness (3 thicknesses 
each for 5CB and 8CB) show no pronounced dependence 
of the undulation wavelength on film thickness, as shown 
in Fig. |4|. In particular, no indication of a quadratic de- 
pendence of the wavelength on film thickness could be 
detected, which would be expected of a spindoal dewet- 
ting scenario. 

It is possible to explain the pronounced effect of the 
heating rate upon the wavelength of the undulation with 
a few assumptions on the nucleation sites. First we note 
that the undulations appear only in a certain temper- 
ature window, Ti < T < T} ll which is more extended 
and shifted to lower temperatures in thin films than in 
thicker ones. For thick films Th will approach the bulk 
N-I phase transition temperature in accordance to earlier 
studies fig] . Each nucleation site present in the film (e. g. 
the domain borders seen clearly in the 8CB samples) is 
expected to nucleate an isotropic domain, and concomi- 
tantly a modulation of the local film thickness, when a 
certain temperature T s within the interval AT = Th — Ti 
is reached. This modulation, or domain then grows with 
a certain velocity v, and thus further nucleation within 
the domain is precluded, effectively reducing the number 
of nucleation sites to become active at a higher temper- 
ature. The total number density of isotropic domains D 
that will develop in a film, and hence the dominant lateral 
scale of the undulation A w l/y/T), thus depend on the 
heating rate. If n(T s ) is the number density of nucleation 
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sites which nucleate a domain when the temperature T s 
is reached, D, will be given by 



dD 
~~dt 



= an{T(t))e- A W 



(1) 



where a is the heating rate, and A(t) is the total area 
fraction of the domains if they are assumed circular. The 
exponential takes care of the mutual overlap of these 
(circular) model domains. For simplicity, assume that 
n(T) = hq = const, for T; < T < Th, and otherwise. 
If we furthermore assume the growth rate dr/dt of the 
domains to be constant, we have 



cPA 
dt 2 



2ttv 2 D 



and thus, after some rearrangements, 
2ttv 2 D(D') 2 + D'D'" — (D") 2 







(2) 



(3) 



For investigating the effect of a temperature ramp, T = 
Ti + at, this must be solved with the initial conditions 
-D(O) = and D'(0) — anQ. The number density of 
domains formed in this way cannot be expressed analyt- 
ically, but is well approximated by 



D„ 



N 



3(a/a ) 2 / 3 
1 + 3(a/a ) 2 / 3 



(4) 



which gives the correct scaling for a — > as well as 
a — > oo. N = no AT is the total number of nucleation 
sites and ao 

does 



V2JWAT2. 
As shown in Fig. |[ the fit function in Eq. 
reproduce the heat rate dependence of the undulation 
wavelength rather well. The scatter in the calculated 
values for N and v is too large to confirm any system- 
atic variation with film thickness, as it is expected for a 
nucleation process triggered by inhomogeneities. A sys- 
tematic dependence on film thickness, as expected for 
spinodal dewetting, is not found. The different evolution 
of the nematic domains in 5CB (created upon spincast- 
ing the films) and 8CB (domains grow at first heating 
cycle), however, could lead to specific differences in the 
two types of samples. Further measurements are neces- 
sary to elucidate this point. 
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FIG. 4. Dependence of undulation wavelength on heating 
rate and film thickness h of 5CB, 51 nm (A), 87 nm (•) and 
139 nm (□). Data for d=87 nm are fitted according to Eq. ^ 
and are presented seperately on a log-log scale in the inset. 

It is finally of interest to investigate the temperature 
range AT, in which the modulation pattern persists. Val- 
ues for AT (based on observation of the temperatures 
at which the undulation started to appear respectively 
ceased to be visible) were found to be independent of the 
heating rate as shown in FigJ| (a) but decrease with 
increasing film thickness, Fig. 0(b). 
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FIG. 5. (a) Temperature interval AT = T h - T t inside 
which the undulations are visible plotted against heating rate 
for three different film thicknesses h, 51 nm (A), 87 nm (•) 
and 139 nm (□). (b) T - T NI (A) and T h ~ T m (o) plotted 
against 1/h 2 . T; and Tj, are mean values over all heating rates 
respectively. The straight lines represent data from Jl7| . Data 
shown are for 5CB films, 8CB films show the same qualitative 
behavior. 

The mechanism leading to different thicknesses of the 
nematic and isotropic parts of the film during the insta- 
bility is still not clear. Initially, it was proposed that 
a Marangoni flow could drive the instability fll8| , since 
the surface tensions of liquid crystalline compounds have 
been shown to exhibit a sharpincrease in the vicinity of 
the N-I phase transition |]l9|-pi|]. Once the phase transi- 
tion is completed, the gradient in surface tension would 
vanish, resulting in a flat isotropic film, explaining the 
transient nature of the instability. However, very recent 
observations of nCB films kept at constant temperatures 
in a temperature range near the bulk phase transition 
temperature have found a nematic-isotropic phase co- 
existence with both phases taking a distinct equilibrium 
thickness |l7j . This effect could very well drive the thick- 
ness modulations in our dynamic experiments. Indeed, 
our plot of AT against hr 2 in Fig. || (b) is in good agree- 
ment with the data found in Jl7j by Effenterre et. al. 
under different experimental circumstances. 

To conclude, we have shown that the undulative in- 
stability observed in thin films of nCB liquid crystals is 
not based on a spontaneous (spinodal) dewetting mecha- 
nism, but develops during the coexistence of the nematic 
and isotropic phase in the films in a temperature win- 
dow close to the N-I phase transition. Isotropic areas 



are nucleated at defects (domain borders) in the nematic 
film and the growth and coalescence of these areas lead 
to the observed undulation patterns. Changing the heat- 
ing rate during the experiments strongly influences the 
undulation patterns since less isotropic areas can be nu- 
cleated at high heating rates before the phase transition 
is complete. First experiments on substrates with ar- 
tificially induced nucleation sites confirm the proposed 
nucleation mechanism. These findings will be covered in 
a forthcoming paper. 
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